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FOREWORD 


This  work  was  performed  for  and  funded  by  the  Independent  Research  Program  at 
Naval  Surface  Warfare  Center  under  IR  task  number  7R01AB313,  Non-Linear  Analysis 
of  Non-Neutral  Plasmas .  This  is  an  investigative  study  undertaken  to  develop  a 
nonlinear  analysis  of  electron  and/or  ion  beams  for  accurate  prediction  of  their 
behavior  as  a  function  of  time.  The  configuration  chosen  for  the  first 
application  of  this  analysis  was  the  TE^  interaction  with  an  axis-rotating 
electron  beam,  in  a  cylindrical  waveguide,  operated  in  the  amplifier  mode.  We 
have  successfully  derived  the  efficiency  of  this  interaction  in  terms  of 
experimentally  controllable  parameters,  and  have  optimized  it  with  respect  to  wave 
frequency. 


Approved  by: 

C.  W.  LARSON,  Head 
Radiation  Division 
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CHAPTER  1 
INTRODUCTION 


GENERAL  BACKGROUND 


The  project,  Nonlinear  Analysis  of  Non-Neutral  Plasmas .  is  ultimately 
concerned  with  the  uses  of  electron  and/or  ion  beams  in  the  generation  of  high 
power  microwave  amd  millimeter  waves  or,  conversely,  the  production  of  high 
brightness  particle  beams.  Military  applications  of  high  power  microwaves  include 
microwave  jamming,  radar  tracking,  and  low  altitude  satellite  neutralization. 

High  brightness  particle  beams  are  under  investigation  as  potential  weapons  for 
application  as  part  of  the  Strategic  Defense  Initiative,  and  may  also  be  used  for 
particle  beam  guidance  systems. 

The  very  high  power/brightness  requirements  of  these  devices  necessitates  the 
development  of  a  new  generation  of  highly  efficient  electron/ion  beam  sources, 
beyond  the  level  of  those  which  are  currently  available.  In  the  past,  most 
attempts  to  build  such  devices  have  relied  on  empirical  data  only,  a  method  which 
has  often  been  costly,  and  which  certainly  carries  no  guarantee  of  success.  It  is 
strongly  felt  within  the  particle  beam  community  that  the  development  of  the  next 
generation  of  electron/ion  beam  devices,  with  their  more  stringent  requirements 
for  power,  gain,  and  frequency,  will  only  be  possible  with  a  much  more  complete 
understanding  of  the  interaction  mechanisms  than  has  been  achieved  to  date. 

These  mechanisms,  involving  the  interaction  of  the  beam  and  the  rf  field, 
differ  from  problem  to  problem,  but  they  have  a  point  in  common:  that  they  are 
usually  intrinsically  nonlinear.  An  examination  of  Figure  3-3  of  this  report 
illustrates  the  point.  In  a  description  of  the  growth  of  the  wave  field  amplitude 
in  a  microwave  generation  problem,  there  is  a  region  of  linear  growth  marked  on 
the  graph,  but  there  is  also  a  region  of  much  larger  nonlinear  growth.  Valuable 
information  about  the  saturation  point  of  the  interaction,  which  fixes  the  maximum 
efficiency  of  the  device,  is  only  accessible  by  studying  the  evolution  of  the  wave 
field  amplitude  through  the  nonlinear  regime. 

Many  examples  such  as  this  one  have  motivated  the  nonlinear  analysis  of 
problems  involving  beam-rf  field  interactions.  There  are  two  approaches  to 
nonlinear  analysis.  One  is  to  perform  a  computer  simulation  of  the  problem.  This 
involves  combining  given  initial  conditions  with  the  electron/ion  equations  of 
motion  and  then  solving  these  equations  in  a  recursive  fashion.  There  are  two 
major  obstacles.  One  is  that  the  number  of  particles  that  even  the  best  modern 
computers  can  handle  falls  far  short  of  the  number  of  particles  in  even  a  low 
current  beam.  Tha  other  is  that  the  evolution  of  the  system  becomes  a  black  box 
with  all  its  basic  physics  hidden.  Optimization  of  results  is  reduced  to  trial- 
and-error  manipulation  of  the  initial  conditions. 
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The  second  approach  to  nonlinear  analysis  is  to  proceed  analytically  as  far 
as  possible,  and  then  to  introduce  approximations  into  the  solution  of  the 
equations  of  motion  in  such  a  way  that  the  important  features  of  the  beam-rf  field 
interaction  are  unaffected.  This  is  the  approach  which  we  have  followed.  Our 
goal  has  been  to  develop  a  nonlinear  analysis  of  electron  and/or  ion  beams,  for 
accurate  prediction  of  their  behavior;  most  importantly,  their  behavior  as  a 
function  of  time.  Our  approach  involves  the  approximate  solution  of  the 
electromagnetic  field  equations,  retaining  nonlinear  features  as  they  pertain  to 
the  time  development  of  electromagnetic  field  quantities  and  beam  parameters. 

This  approach  has  proved  very  successful  in  the  analysis  of  the  first 
configuration  chosen  for  examination,  the  cylindrical  waveguide  TE^  interaction 
with  an  axis-rotating  electron  beam,  operating  in  the  amplifier  mode. 


CONFIGURATION:  TEn  INTERACTION  WITH  AXIS -ROTATING  BEAM 

The  TE-q  mode  excitation  is  of  great  practical  interest  in  microwave 
generation  since  it  is  the  dominant  mode  in  a  cylindrical  circuit.  The  optimum 
beam  configuration  for  exciting  the  TE^  mode  is  known  to  be  an  axis-rotating 
beam1,  which  can  be  produced  by  injection  of  a  hollow  beam  through  a  cusp 
magnetic  field.  In  preliminary  experiments  powerful  radiation  with  multikilowatt 
output  at  the  fundamental  electron  cyclotron  frequency  was  observed  in  the  NSWC 
Cusptron  device.  This  prompted  a  detailed  theoretical  investigation  of  the 
efficiency  of  this  configuration. 


APPROACH 


To  theoretically  study  the  efficiency  of  TE-q  mode  generation  by  an  axis- 
rotating  beam  via  the  negative  mass  instability,  a  nonlinear  analysis  is 
developed.  Other  configurations  analyzed  by  nonlinear  analysis  include  the 
gyrotron  in  a  cylindrical  waveguide i  and  in  a  parallel  plate  waveguide4;  and  an 
axis -rotating  beam  in  a  multivane  structure^.  We  assume  that  phase  trapping  of 
the  electrons  in  the  wave  fields  is  the  dominant  saturation  mechanism.  The 
investigation  of  the  problem  starts  with  solving  the  nonlinear  coupled  equations 
of  single  electron  motion  in  the  radiation  fields.  Next,  the  response  of  the 
radiation  fields  to  the  current  induced  by  the  electron  bunching  is  analyzed  using 
conservation  of  energy.  Then  the  results  are  averaged  over  the  initial  random 
phase  distribution  of  the  electron's  transverse  velocity  to  obtain  the  collective 
effect  on  the  wave  fields  by  the  electrons.  The  result  of  the  nonlinear  analysis 
performed  in  this  study  is  the  derivation  of  a  single  nonlinear  third-order 
differential  equation  for  the  wave  field  amplitude  as  a  function  of  time.  From 
knowledge  of  the  temporal  evolution  of  the  wave  field,  the  efficiency  of  the 
interaction  may  be  directly  calculated. 
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REPORT  OVERVIEW 


In  Chapter  2  we  develop  the  nonlinear  analysis  necessary  to  study  the  wave 
field  amplitude  as  a  function  of  time.  We  derive  a  third-order  differential 
equation  which  may  be  solved  numerically  to  obtain  the  temporal  evolution  of  the 
wave  field  over  its  entire  range,  through  the  linear,  nonlinear,  and  saturation 
regimes.  In  Chapter  3  the  numerical  solution  of  the  full  equation  is  carried  out 
and  the  efficiency  of  the  interaction  is  studied  parametrically  as  a  function  of 
the  beam  energy  and  beam  current,  the  external  magnetic  field,  and  the  waveguide 
radius.  The  efficiency  is  also  optimized  as  a  function  of  wave  frequency.  Our 
conclusions  are  presented  in  Chapter  4. 
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CHAPTER  2 

NONLINEAR  ANALYSIS  OF  TEn 
INTERACTION  WITH  AXIS -ROTATING  ELECTRON  BEAM 


FIELDS  AND  GEOMETRY 


The  geometry  of  the  problem  is  as  shown  in  Figure  2-1.  The  equilibrium 
configuration  consists  of  a  non-neutral  electron  layer  that  is  infinite  in  extent 
and  aligned  parallel  to  an  applied  magnetic  field  B  z.  The  electron  layer  is 
accomplished  by  passing  a  hollow  beam  through  an  ideal  cusp  magnetic  field.  The 
investigation  of  the  problem  starts  with  solving  the  nonlinear  coupled  equations 
of  single  electron  motion  in  the  electromagnetic  fields.  The  motion  of  an 
electron  in  the  external  magnetic  field  and  the  waveguide  mode  fields  is  governed 
by  the  set  of  coupled  nonlinear  equations 


dr 


dt  ym  ’ 
o 


— e [ E  +  ±  v  x  (B  +  Bqz)]  , 


dt 
G7  e 


dt  2 

m  c 
o 


E  *v  , 


(2-1) 

(2-2) 

(2-3) 


where  7  -  (l  +  P^/(m^c^)J^  ,  and  P 


7m  v  . 
o 


The  wave  fields  E(x,t)  and  B(x,t)  of  the  TEmn  mode  inside  a  circular  cylinder 
of  radius  Ry  are  given  by, 


mJ  (7  r) 

E  --E  (t)  — — — — —  cos(m<£  +  k  z  -  wt)  , 
r  o  7  r  ^  z 


mn 


E ,  -  E  (t)  J '  (7  r)  sin(m<£  +  k  z  -  wt)  , 
<f>  o  m  mn  z 


E  -  0  , 
z 


k  c 

B  — —  E  , 
r  w<f> 
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7  c 

B  =  -mn  E  (t)  J  (7  r)  cos(m4>  +  k  z  -  wt)  ,  (2-4) 

z  w  o  m  mn  z 

where  7^  is  the  n*"*1  root  of  the  equation:  J^(7mnR^)  =0.  We  assume  that  the 
electromagnetic  fields  are  those  of  a  vacuum  waveguide,  but  with  amplitudes  that 
vary  with  time.  The  implications  of  this  assumption  are  discussed  in  Chapter  4. 

Rather  than  approach  this  problem  as  a  computer  simulation  we  will  try  to 
identify  the  source  of  the  greatest  contribution  to  the  wave  field  growth.  The 
exchange  of  energy  between  the  wave  field  and  the  beam  will  be  greatest  when  the 
wave  frequency,  w,  is  near  the  Doppler  shifted  electron  cyclotron  frequency, 
mfl  +  k  v  ,  where  fl  —  0  /y,  and  fl  =  eB  /m  c.  In  this  paper  we  examine  the 
fundamental,  m  -  1,  case.  We  note  that  tRe  quantity 


Au(t)  =  w  -  {kzvz(t)  -  fl(t)} 


(2-5) 


essentially  defines  the  degree  of  resonance  present.  In  this  approach,  an  average 
is  taken  over  the  initial  random  phase  in  such  a  way  that  quantities  related  to 
the  fast  cyclotron  motion  are  averaged  out. 

A  self-consistent  trajectory  is  employed  to  describe  the  location  of  a  single 
electron  within  the  rotating  electron  layer, 


r(t)  -  RL(t)  , 

(2-6) 

*(t)  =  4>Q  +  *(t)  +  /J  n(t’)dt'. 

(2-7) 

RL(t)  is  the  time -dependent  Larmor  radius  of  the  rotating  electron.  The  time- 
dependent  phase  4>( t)  is  composed  of  three  parts:  the  initial  phase  <f>  ,  the  usual 

term  due  to  the  electron  cyclotron  motion  in  the  external  B  z  field,  and  a  slowly 
varying  phase  ip(t)  due  to  the  interaction  of  the  wave  field  with  the  electron.  By 
substituting  Equations  (2-1),  (2-4),  (2-5),  and  (2-6)  into  Equations  (2-2)  and  (2- 
3),  a  total  of  four  equations  are  derived  for  the2four  slowly  varying  quantities: 
the  total  relativistic  factor  y(t)  (energy  =  7m  c  ),  the  axial  velocity  v  (t),  the 
slowly  varying  phase  'f(t) ,  and  the  generalized  Rarmor  radius  a(t)  -  y^^R^ft) . 
Explicitly,  these  equations  are 
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da(t) 

dt 

d$(t) 

dt 


where  f (a)  -  aJ'(a)/y  ,  c-  =  711c/(B  w) ,  and  c,  -  c„(w  -  k  v  ) 
l  l  11  o  l  z  z  z 

The  quantity  $>(t)  is  the  relative  ph^se  $>(t)  -  <£(t)  +  k  z(t)  -  wt  ,  which 
may  be  expressed  as  $(t)  =  (<£  +  V»( t ) )  -  /  Aw(t')dt'  where  Xw(t)  is  given  by 

Equation  (2-5).  °  U 


-E  (t)c1Ji (a)sin  $  , 
oil 


(2-10) 


-  Aw ( t ) 


(2-11) 


*?v>  • 


cos  $>  , 


GRAZING  CONDITION 


The  efficiency  is  expected  to  be  highest  at  or  near  grazing,  therefore  we 
limit  our  attention  to  the  case  where  the  grazing  condition,  v  =  k  c  /w  ,  is 
valid.  Upon  application  of  the  grazing  condition,  Equation  (2-8)  for  dvz(t)/dt 
equals  zero;  thus  the  grazing  condition  implies  that  v2(t)  -  vZQ  at  all  times. 

From  Equation  (2-11),  Aw(t)  may  be  explicitly  evaluated.  This  is  easily  shown  to 
be 

,2  2 

k  c  7-7 

Aw(t)  -  Aw  +  fw  -  — - ]  f - -}  ,  (2-12) 

O'-OW-'V-yJ 


where  Aw  —  w  - 


and  w  =0/7  +  k  v  ,  although  the  grazing  condition  must  be 


O  O  O  D  O  2  Z  Q  0  00 

used  to  obtain  the  result  in  this  form.  Integrating  Equation  (2-8)  to  obtain 


7  -  7, 


Equation  (2-12)  becomes 


Aw ( t )  —  Aw 


,2  2 
k  c 
z 

w 


)  Jq  dt'E0(t')  wc2  [^)2f(t')  Sin  *(t')}  .  (2-13) 


CONSERVATION  OF  ENERGY 


At  this  point,  Equations  (2-7) -(2-10)  define  the  effect  of  the  wave  field  on 
a  single  electron.  To  find  the  reciprocal  effect  of  the  ensemble  of  beam 
electrons  on  the  wave  field  an  additional  field  equation  is  required. 

Conservation  of  energy  (Poynting  equation),  with  the  fields  in  Equation  (2-4)  is 
adequate  for  this  purpose  as  Equations  (2-7) -(2-10)  involve  only  the  amplitude  of 
the  wave  field.  It  is  straightforward  to  show  that  conservation  of  energy  implies 
that 
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n  m  c 
o  o 


2/d7(t)  V 

dt  / 

\  / 


+ 


E  (t) 
o 


diyt) 

dt 


=  0  . 


(2-14) 


where  nQ  is  the  electron  beam  density  averaged  over  the  cross-section  of  the 
waveguide,  <  >  stands  for  the  average  over  the  initial  random  phase  distribution 

of  the  electron  beam,  and  e  is  given  by 


-  t1  -  t-H  • 

711RV 


(2-15) 


Upon  substitution  of  Equation  (2-7)  for  dy(t)/dt  into  Equation  (2-14)  one  finds 
that 


dE  (t) 

<fsin($)>  =  — 


P  dt 


(2-16) 


where  P  =  47rn  m  Q  c/(f7,  ,B  ) 
o  o  o  '11  o 


PHASE  AVERAGE 


We  next  determine  the  phase  averaged  quantity  <fsin  $  >  in  terms  of  EQ(t)  and 
its  derivatives  using  the  two  remaining  equations,  Equations  (2-10)  and  (2-11). 
Details  of  the  calculation  are  provided  in  Appendix  A.  Three  approximations  are 
made.  First,  it  is  assumed  that  the  effect  of  beam  energy  depletion  can  be 
neglected.  This  is  consistent  with  our  assumption  that  phase  trapping  is  the 
dominant  saturation  mechanism.  Second,  it  is  assumed  that  wave  frequency  is  on 
the  order  of  (but  not  equal  to)  the  initial  electron  cyclotron  frequency.  The 
third  approximation  is  that  only  linear  and  nonlinear  terms  of  the  lowest  order 
are  retained  in  the  performance  of  the  phase  averaging.  Details  of  the  phase 
averaging  procedure  are  provided  in  Appendix  B  and  discussed  in  Chapter  4. 


NONLINEAR  DIFFERENTIAL  EQUATION  AND  EFFICIENCY 


The  final  result  is  now  derived.  Using  Equation  (A- 7)  for  <fsin  $  >  and 
Equation  (A-8)  for  <fcos  $>  >  ,  together  with  Equation  (2-16),  the  following 
nonlinear  third-order  equation  for  the  wave  field  amplitude  EQ(t)  as  a  function  of 
time,  t,  results, 
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where 
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and 


Is  =  -/£  d£'X(£')sinj|,[l  -  2g(XZ(|)  -  XZ(0))]d?  , 


2  - 


Ic  =  d^X(r)cos/|,[l  -  2g(X2(0  -  XZ(0))]d? 


(2-21) 


The  efficiency  of  the  interaction  may  now  be  calculated  according  to^’^ 


2f  X2  -  X2  ) 

1  -  — T-5 - ’  (2-22) 

°1 

where  Xs  is  the  value  of  the  normalized  electric  field  amplitude  at  the  point  of 
saturation.  In  Chapter  3,  Equations  (2-19)  and  (2-22)  will  be  solved  numerically 
to  determine  the  optimized  efficiency  as  a  function  of  various  parameters:  beam 
energy  and  beam  current,  external  magnetic  field,  and  waveguide  radius. 
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FIGURE  2-1.  EQUILIBRIUM  CONFIGURATION  AND  COORDINATE  SYSTEM. 
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CHAPTER  3 
NUMERICAL  RESULTS 


The  TE^  interaction  of  fast  axis-rotating  electron  beams  in  a  cylindrical 
waveguide  may  now  be  analyzed  using  Equation  (2-19)  to  determine  the  evolution  of 
the  wave  field  through  the  linear  and  nonlinear  regimes  up  to  the  point  of 
saturation,  and  Equation  (2-22)  to  evaluate  the  corresponding  efficiency. 

We  choose  the  initial  conditions  from  among  the  experimentally  relevant 
parameters.  These  are  the  beam  energy  E^,  the  beam  current  I,  the  external 
magnetic  field  BQ,  and  the  waveguide  radius  Rw. 

It  is  important  to  solve  for  the  efficiencies  within  the  limits  of  the 
assumptions  used  to  produce  Equation  (2-19).  We  examine  the  consequences  of  our 
assumption  that  the  grazing  condition  is  always  met  before  proceeding  to  the 
efficiency  calculation. 


CONSEQUENCES  OF  GRAZING  CONDITION 


The  grazing  condition  is  satisfied  at  the  point  where  the  beam  dispersion 
curve  tangentially  grazes  the  wave  dispersion  curve  , 


with 


k  v 
z  zo 


r  2  2.  2 

to  -  I  v  +  c  k 
w  t  co  z 


to 

CO 


5.523  x  1010 


i/2- 

/  ^  • 


(3-1) 

(3-2) 


for  the  TE^  interaction.  Grazing  occurs  when 
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zo 
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(”2  •  Hr) 

v  co  2  > 


n 


co 


(3-3) 


Writing  this  in  terms  of  Bq  ,  E^,  and  R^  leaves  the  grazing  condition  expressed  in 
terms  of  the  chosen  parameters, 


NSWC  TR  87-340 


2 

a 


R\flo^2 
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(3-4) 


where  a  =  v  /v  This  equation  is  only  true  for  certain  combinations  of  the 
parameters  E^,  B0,  and  Ry.  In  any  solution  of  Equation  (2-19)  for  the  wave  field 
amplitude  X(£)  (the  normalized  EQ(t)),  a  is  first  determined  for  a  given  Ry  and  BQ 
so  that  the  grazing  condition,  Equation  (3-4) ,  is  satisfied. 

Figure  3 "3-  shows  the  behavior  of  a  versus  R^/Ry,  along  with  the  corresponding 
values  for  P  and  p  .  This  behavior  has  implications  for  the  efficiency 
calculation.  Although  a  can  asymptotically  approach  infinity,  it  is  clear  that 
there, will  be  no  increase  in  efficiency  beyond  the  point  where  P  goes  to  zero 
and  p  reaches  its  maximum,  i.e.,  the  point  where  all  the  initial  energy  is 
assumed  to  go  into  rotational  energy.  Therefore,  while  RW  and  B  may  both  be 
independently  varied,  application  of  the  grazing  condition  will  link  them  together 
and  limit  the  effect  of  their  individual  variation  on  the  efficiency  to  a  maximum 
at  the  point  where  p  -  0  .  This  point  is  illustrated  in  Figure  3-2,  a  plot  of 
efficiency  versus  a  for  a  30.0  keV,  3.5  Amp  beam  in  a  460.0  Gauss  external 
magnetic  field. 


OPTIMIZED  EFFICIENCIES 


We  may  now  numerically  solve  Equation  (2-19)  in  a  self-consistent  manner.  A 
30.0  keV,  3.5  Amp  beam  in  a  460.0  Gauss  external  magnetic  field  is  chosen  for 
study.  Additional  studies,  to  be  reported  elsewhere,  indicate  that  changes  in  E^ 
and  I  do  not  greatly  change  the  results  to  be  presented  here.  Variations  in  BQ 
and  Ry  have  been  discussed  in  the  previous  sub-section;  for  a  460.0  Gauss  field, 
an  Ry  -  7.227  cm  results  in  the  best  possible  efficiency. 

A  Runge-Kutta  scheme  was  employed  to  evaluate  X(£)  versus  £  ,  and  for  most 
combinations  of  parameters  a  step  size  of  -  0,0001  was  adequate.  A  single 
efficiency  calculation  involving  2000  steps  takes  approximately  3.953  CP  seconds 
on  a  CDC  CYBER  720  mainframe  computer. 

In  Figure  3-3  we  display  a  typical  result  obtained  by  integrating  Equation 
(2-19).  It  shows  the  evolution  of  the  wave  field  amplitude  for  a  30.0  keV,  3.5 
Amp  beam  in  an  external  magnetic  field  of  460.0  Gauss  with  Ry  -  7.227  cm.  The 
linear,  nonlinear,  and  saturation  regions  are  indicated  on  the  graph.  The  point 
of  saturation  is  clearly  defined  at  X  -  0.0517.  The  corresponding  efficiency  is 
9.11%. 

It  is  desired  to  perform  not  only  a  single  efficiency  calculation,  but  an 
efficiency  optimization  study.  We  choose  to  optimize  the  efficiency  with  respect 
to  wave  frequency  as  that  is  an  experimentally  variable  parameter.  To  determine 
an  appropriate  parameter  range  for  w  we  utilize  earlier  results  from  the  linear 
theory  for  this  configuration'*  The  result  of  this  analysis  is  an  optimized 


VCWiA.'A.AV.V-V.V.V/i’.V.Wirt.V-V 
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efficiency  as  a  function  of  w  . 

In  Figure  3-4,  a  plot  of  efficiency  ry  as  a  function  of  frequency  f^ 

(_f  -  w/2 n  )  is  shown.  The  efficiency  is  clearly  optimized  at  f  *  1.233  GHz,  with 
an  efficiency  of  9.11%. 

The  curve  itself  is  likely  t£  be  artificially  broadened  by  the  necessary  use 

of  an  approximation,  k  “  w  v  /c  ,  in  regions  where  Aw  is  not  small.  Studies 

z .  p .  zp  ,  o 

are  currently  underway  in  which  the  grazing  condition  assumption  itself  is 
removed. 


r 

- 


! 


► 
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Rl/Rw 


FIGURE  3-1.  BEHAVIOR  OF  a ,  AND  0^  AS  FUNCTIONS  OF  RL/Ry 
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NONUNEAR  SATURATION 


FIGURE  3-3 .  EVOLUTION  OF  THE  WAVE  FIELD  AMPLITUDE.  (FOR  A  30.0  keV , 

3.5  AMP  BEAM  IN  A  460.0  GAUSS  EXTERNAL  MAGNETIC  FIELD,  THE 
SATURATION  AMPLITUDE  Xg  IS  0.0517,  CORRESPONDING  TO  A  9.11% 
EFFICIENCY.) 
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CHAPTER  4 

SUMMARY  AND  CONCLUSIONS 


SUMMARY  OF  RESULTS 


We  have  derived  a  single  nonlinear  third-order  differential  equation 
(Equation  (2-19))  to  describe  the  evolution  of  the  wave  field  amplitude  in  the 
TE-j^  interaction  with  a  fast  axis  -  rotating  electron  beam.  From  knowledge  of  the 
evolution  of  the  wave  field  amplitude  the  efficiency  of  the  interaction  may  be 
directly  calculated.  We  assume  phase  trapping  to  be  the  dominant  saturation 
mechanism.  The  derivation  is  self-consistent,  with  the  effect  of  the 
electromagnetic  fields  on  the  electrons  (phase  averaged  over  the  initial  random 
distribution  to  obtain  a  collective  effect)  coupled  to  the  response  of  the  wave 
fields  to  the  current  induced  by  the  electron  bunching  through  conservation  of 
energy . 

Numerical  solution  of  Equation  (2-19)  yields  the  wave  field  amplitude  as  a 
function  of  time  over  its  entire  range,  through  the  linear,  nonlinear  and 
saturation  regimes.  The  efficiency  is  directly  calculated  as  a  function  of  the 
experimental  parameters  beam  energy,  beam  current,  external  magnetic  field,  and 
waveguide  radius.  The  efficiency  is  also  optimized  with  respect  to  wave 
frequency,  with  a  resulting  best  efficiency  of  9.11%  at  a  frequency  of 
f  =  1.233  GHz  . 


APPROXIMATIONS  AND  ASSUMPTIONS 

It  is  interesting  to  consider  the  implications  of  some  of  the  approximations 
and  assumptions  made  throughout  the  analysis.  These  are: 

(a)  .  The  electromagnetic  fields  are  taken  to  be  vacuum  waveguide  fields  with 
time  dependent  amplitudes. 

The  choice  of  vacuum  waveguide  electromagnetic  fields  with  time -dependent 
amplitudes  implies  that  the  electron  beam  itself  is  tenuous;  otherwise  space 
charge  effects  would  considerably  distort  the  electromagentic  fields  from  their 
vacuum  waveguide  forms.  This  is  a  reasonable  assumption  to  make  for  a  30.0  keV, 
3.3  Amp  beam,  but  would  have  to  be  modified  in  a  high  current  situation. 

(b) .  Only  linear  and  nonlinear  terms  of  the  lowest  order  are  retained  ir.  the 
performance  of  the  phase  averaging. 

That  the  phase  does  not  evolve  reciprocally  with  the  growth  of  the  wave  field 
is  a  pronounced  drawback.  We  are  currently  exploring  methods  to  include 
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higher  order  terms  in  our  future  work. 

(c).  Imposition  of  the  grazing  condition. 

For  the  amplifier  study,  the  grazing  condition  is  expected  to  be  valid  at  the 
point  of  greatest  efficiency,  but  not  elsewhere.  The  efficiency  optimization 
study  indicates  that  this  is  indeed  the  case.  More  tantalizing  still  are 
indications,  both  theoretical  and  experimental,  that  slightly  off-grazing 
conditions  may  produce  much  higher  efficiencies.  The  next  extension  of  this  work 
is  to  remove  the  grazing  condition  requirement.  An  additional  advantage  is  that 
an  oscillator  study  will  then  be  possible,  for  which  experimental  data”  is 
available  for  direct  comparison. 


CONCLUSIONS 


The  accomplishments  of  FY87  are  three.  We  have  successfully  developed  the 
complete  nonlinear  analysis  of  the  TEn^  interaction  with  an  axis  -  rotating 
electron  beam,  in  a  cylindrical  waveguide,  operated  in  the  amplifier  mode.  We 
have  performed,  for  the  first  time,  a  direct  efficiency  calculation  for  theTE^ 
mode  interaction.  We  have,  also  for  the  first  time,  carried  out  an  optimization 
with  respect  to  wave  frequency,  an  experimentally  controllable  parameter. 

It  should  be  remarked  that  the  present  method  offers  significant  advantages 
in  the  numerical  work  involved  as  opposed  to  an  undiluted  particle  simulation,  by 
averaging  out  the  fast  time  scale  part  of  the  Lorentz  force.  It  is  also 
anticipated  that  analytical  studies  of  Equation  (2-19),  currently  underway,  will 
yield  insight  into  the  nature  of  the  physical  processes  taking  place. 


While  the  project  was  particularly  directed  toward  the  analysis  of  the  TEn^ 
interaction  in  FY87,  the  formulation  up  to  Equation  (2-5)  is  general.  The  results 
shown  in  this  report  are  easily  extended  to  the  TEn^  case.  Studies  of  the  TE21 
and  the  TEj^  interactions  are  in  progress  at  the  present  time. 
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APPENDIX  A 

EVALUATION  OF  <fsin  $  > 


We  determine  the  phase  averaged  quantity  <fsin  $  >  in  Equation  (2-16)  in 
terms  of  E  (t)  and  its  derivatives  using  Equations  (2-10)  and  (2-11).  To  include 
adequate  information  provided  by  the  two  coupled  independent  first-order 
differential  equations  a  second-order  differential  equation  is  required.  We 
therefore  examine  d^fsin  $J/dt  and  d^fsin  $  J/dt  .  Using  Equations  (2-10)  and 
(2-11)  to  evaluate  4,  4,  a,  and  a  , 


d^fsin  <f>  ) 
dt 


-E  c 
o  1 


2.  ,  1  d(QJ0  2.  ^ 
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(A- 1 ) 
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dt 

-  fcos  $  +  Eqc^3Au  [f'Jf  ‘  q  da  }sin  $  cos  $ 

+  E2c2{  (f"Jf2  +  f'JfJf]sin3$ 

+  f(Ji(  k  6  k  (oji0  )  -  2S  £(“Ji>]2)sln  * cos2t  }  ■  <A-2> 

£"  -  2c3JJ/i3  +  <o/l){jJ  •  3o2(no/-»n)2(jy14)}  , 

The  dots  •  and  primes  '  refer  to  differentiation 
with  respect  to  time  t  and  a  respectively. 

Two  approximations  are  made.  First,  it  is  assumed  that  the  effect  of  beam 
energy  depletion  can  be  neglected.  This  is  consistent  with  the  assumption  that 
phase  trapping  is  the  dominant  saturation  mechanism.  The  consequences  of  assuming 
that  | A-y ( t )  |  «  7  are:  that  y(t)  -  y  and  a(t)  =  a  .  Second,  it  is  assumed 
that  the  wave  frequency  is  on  the  order  of  (but  not  equal  to)  the  initial  electron 
cyclotron  frequency:  w  '  w  .  With  these  two  approximations,  Equations  (2-5),  (A- 
1),  and  (A-2)  become, 


where  f'  -  aJ'^/y  +  c-Jf/y3, 
and  c3  -  l/{l  -  (vzo/c)2}  . 
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where  c  =  f - - — ]2  c  (w2  -  k2c2)  ,  and  f  =  f(a  ),  J,  -  J, (a  ) ,  etc. 
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The  phase  averages  over  the  initial  random  phase  distribution  of  the 
electron's  transverse  velocity  are  taken  to  obtain  the  collective  effect  of  the 
wave  field  on  the  electrons.  In  Appendix  B,  explicit  expressions  for  the  phase 
averages  of  the  combinations  of  trigonometric  functions  in  Equations  (A-3),  (A-4), 


A-  2 
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and  (A- 5)  are  derived.  A  third  approximation  is  made:  only  the  linear  and 
nonlinear  terms  of  the  lowest  order  are  retained  in  the  performance  of  the  phase 
averaging.  Substituting  these  approximate  phase  averages  into  Equations  ( A -  3 ) , 
(A-4) ,  (A-5)  and  carrying  out  the  necessary  integrations  one  finds  that,  together 
with  Equation  (2-16), 
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4  W07llcJ  31  VllcJ  °  10 


4J 


lo  da 


)  >■ 


I  -  dt'E  (t')  cos  <A$(t  -  t')>  , 
c  J  0  o 


I  -  dt'E  (t')  sin  <A*(t  -  t')>  , 
s  J  U  o 


(A-  9 ) 


with  <A$(t  -  t')> 


» <Aw( t" )>dt" . 

A-  3 
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APPENDIX  B 


EVALUATION  OF  PHASE  AVERAGES 


By  using  the  half  and  double  angle  formulas  we  obtain  the  following  exact  and 
approximate  expressions  for  the  phase  averages: 


)  1  /•  \ 
"$  >  =  <  —  |^1  -  cos  2$  J> 


2  1 

<cos  $  >  =  <  —  (1  +  cos  2$  )  > 


<sin  $>  cos  <1>  >  -  sin  2$  > 


<sin^<$  >  ”  <4  (^sin  $  +  sin  3$  ^  <sin  4>  >  , 


2  r  2  \ 

<sin  3>  cos  $  >  =  <sin  $  II  -  sin  $  J> 


y  <sin  $  >  . 
4 


(B-l) 


Higher  order  terms  involving  <cos(2«I>)>,  <sin(2<S>)>,  and  <sin(34>)>  have  been 
neglected. 


Angle  addition  and  subtraction  laws  are  used  to  evaluate  the  following 
expressions 


<sin  $>  sin  $'  >  =  |cos($-$' )  -  cos ($+$'  )^>  *<  ^  cos<A$(t  -  t')>  , 

<cos  $  sin  $'  >  =  <^  ^-sin($-$')  +  sin($+<I>' ) j->  =  -^  sin<A$(t  -  t')>  , 

<sin  $  cos  $  sin  $'  >  =  ^  {<sin  $  (cos  $  sin  $'  )>  +  <cos  <E>  (sin  $  sin  $')>^  , 


<sin  $  >sin<A$(t  -  t')>  +  ^  <cos  $  >cos<A3>(t  -  t')>  .  (B-2) 


Higher  order  terms  involving  <cos(3>  +  $')>  and  <sin($>  +  $'  )>  have  been  neglected 


for  the  following  reason.  Based  on  Equation  (2-11),  $  may  be  generally  expanded- 
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where  <$>  -  <$ 

o 


2 

$  -  <4’>  +  C^sin(<$>  +  0^)  +  C^sin  (<<t>>  +  ©2)  +  . 
-  Jq  <Aw( t ' )>dt ' .  Therefore 

rt 

$  -  $>'  --J  <Au(t")dt"  +  higher  order  terms 
$  +  $'  -  higher  order  terms  , 


thus  we  approximate 

<$-$'>=  <A$(t  -  t')>  , <Aw( t" )>dt" . 


and  neglect  terms  involving  cos<  $  +  $'  >  and  sin<  4>  +  >. 
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